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ABSTRACT 



We studied both components of a slightly overlooked visual binary HR 1847 spectroscopically to determine its basic physical and 
orbital parameters. Basic stellar parameters were determined by comparing synthetic spectra to the observed echelle spectra, which 
cover both the optical and near-lR regions. New observations of this system used the Ondfejov and Rozhen 2-m telescopes and their 
coude spectrographs. Radial velocities from individual spectra were measured and then analysed with the code FOTEL to determine 
orbital parameters. The spectroscopic orbit of HR 1847 A is presented for the first time. It is a single-lined spectroscopic binary with 
a B-type primary, a period of 719.79 days, and a highly eccentric orbit with e = 0.7. We confirmed that HR 1847B is a Be star. Its Ho- 
emission significantly decreased from 2003 to 2008. Both components have a spectral type B7~8 and luminosity class IV- V. 
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1. Introduction 

; HR1847 (HD 36408, BD+16°794, HIP 25950, ADS 4131) is a 
bright visual binary consisting of two B type stars. Throughout 
this paper, the component at a'(J2000) = 05 32 14.14, 5(J2000) = 
+ 170329. 3 will be denoted as HR 1 847 A (or component A) and 
. the one at q'(J2000) = 05 32 14.56, 5(J2000) = +17 03 21.8 as 
■ HR 1847B (or component B). The first observations of this bi- 
nary, which were obtained on 31 December 1782, were reported 
. by William Herschel (see Herschel, 1785, star III.93). He ob- 
' served this binary again on 2 1 January 1 800 (see lHerschelll 1 82 U 
' star 124). 

Athough th is visual binary consists of two bright stars and 
, iPlaskettI (Il919l) recommended further observations, it was not 
I studied much during the last centur y. Since HR 1847 is both 
[ an X-ray source (Berghofer, Schmitt, & Cassinellill 19961) and an 
IRAS source (Cote & Waters, 1987), and considering its bright- 
ness, lack of observations, and sometimes confusing catalogue 
information, it has become quite an interesting object for more 
! detailed study. 



2.1. Multiplicity 



2.1.1. Angular separation 



Angular separation p and relative position 6 of the components 
of HR 1 847 have alrea dy been measured by Friedrich Georg 
WiUielm IStruvd (1 18371 star No. 730) as an average of 4 mea- 
surements obtained in 1829 and 1832 (p = 9'.'81,6I = 141°82). 
Later m easurements resulted in p = 9753,0 - 140?2 (|CannonL 
IT9T2bl) . o = 9'.'7. 6> ^ 140°9 (iRussell & Moord [T929I). o ^ 
9:65,9 = 140°85 dPerrvman et al.L 1 19971 from lllPPARCOSl. 
p = 9'.'65 + 0'.'006, = 140°5 ('Fabricius et al.', '2002' from the 
Tycho catalogue), p = 9'.'61, = 140?8 (Mason et al, 2004, from 
speckle interfero met ry). The angular separation, as reported by 
iLindroosI (119831) and lAbt & Cardonal(il984) is 9'.'6. 



2.1.2. Radial velocity variations 



2. Summary of known properties of HR 1847 

In this section we summarize the information currently available 
for both components of HR 1847, which is scattered throughout 
the astronomical literature, online catalogues, and databases. 



Send offprint requests to: J. Kubat, 
e-mail: kubatdsunstel . asu . cas . cz 

* Based on observations with the Ondrejov and Rozhen 2-m tele- 
scopes. 



The first radial velocity measurements for both components were 
rep orted by [Ha skett & Youn^ (Il919l) and they showed variabil- 
ity. iPlaskett et al.i ([1920) suggested that both components are 
spectroscopic binaries. [Frost et al. ( 1926) designated HR 1847A 
(No. 126 in his Table II) as a spectroscopic binary with an un- 
known orbit using only one spectrum secured at the Dominion 
Astrophysical Observatory (DAO) in 1919 by Otto Struve. 
Component A was also reported to be a spectroscopic binary 
by Gahm. Ahlin. & Lind roos (1983|) andJAbt. Levato. & Grossj 
(l2002h . However, no attempt has been made to determine the 
orbital parameters until now. 
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Table 1. HIPPARCOS and UBVRIJHK photometry of 

HR1847. 





HR1847A 


HR1847B 


source 




band 


mag 


mag 




Hn 


6.088 ± 0.007 


6.455 ±0.010 


HIPPARCOS 




Br 


6.036 ± 0.008 


6.468 ± 0.004 


HIPPARCOS 




Vr 


6.057 ± 0.009 


6.458 ± 0.004 


HIPPARCOS 




U 


5.823 ± 0.030 


6.337 ± 0.030 


EXPORT 




B 


6.103 ±0.037 


6.527 ± 0.037 


EXPORT 






6.070 ± 0.020 


6.520 ± 0.020 


Maidanak 




V 


6.090 ±0.100 


6.490 ± 0.077 


EXPORT 






6.090 ± 0.020 


6.510 ±0.020 


Maidanak 




R 


6.037 ±0.100 


6.410 ±0.077 


EXPORT 






6.050 ± 0.020 


6.460 ± 0.020 


Maidanak 




I 


6.050 ±0.150 


6.437 ±0.133 


EXPORT 




J 


5.964 ±0.019 


6.280 ± 0.024 


2MASS 




H 


6.022 ±0.021 


6.275 ±0.021 


2MASS 




K 


6.016 ± 0.027 


6.282 ±0.016 


2MASS 




HIPPAP 


.COS - Perrvma 


n et al. ( 1997), EXPORT - Oudmaijer et al.' 


(2001), 


Vlaidanak -Shats 


kii(1998),2MA. 


^S -.Skrutskie et al. (2006!). 



2.2. Distance 

The HIPPARCOS par allax of HR1847 is 2.92 + 1.57 mas 
dPerryman et al.L |1997|) . which corresponds to a distance of 
342 + 184pc, but the large error reduces the reliability of this 
distance es timate. However, the new edition of the HIPPARCOS 
catalogue dvan Leeuwenl 120071) gives a parallax value of 1.09 + 
2.04 mas, which is unable to provide any information about the 
distance owing to the very uncertain parallax. Individual parallax 
measurements of the A and B components are not available. 



2.4. Polarization 

lOudmaiier et al.l (1200 lb found the value of V band polarization 
Pv = 0.0063 + 0.0004 for the A component, and Py = 0.0063 + 
0.0002 for the B component. 

2.5. Spectral types and luminosity classes 

The early class ification of bo th comp onents w as 'e^regie alba' 
(very white) bv lStruvel(ll837h . Later, (CannonL[l912b) classified 
both stars as B9. 

Using his own observations in the period from April 1954 
to March 1955 at the Yerkes Observatory, Osawa (1959) deter- 
mined the spectral type of component A as B7IV and of compo- 
nent B as B8IV. From his own observations with the Perkins 
telescope, ISlettebakI (Il963h determined spectra l types B7III 
for component A and B7IV for component B. ICowlev et al.l 
(1969) determined spectral types B8III and B8V for compo- 
nents A and B, respecti vely, usmg sp ectra from MacDonald 
and Yerkes Observatories. lLevatol (1197 5) observed both compo- 
nents at Cerro Tololo Inter-American Observatory and obtained 
spectral type B7III for component A and B8IV for B. Finally, 
[Mora et al. (2001) changed the Bright Star Catalogue classifica- 
tion of HR 1847 A from B7IIIe to B5V. 



2.5.1. Emission 

[Andrews! (1 1968h found that component B shows Ha emission us- 
in g Ha photornetry, \y hich was later spectroscopically confirmed 
bv'Wack erling|(ll970h . The emission sign, which appeared at the 
spectral type of component A in both the Bright Star Catalogue 
and the Simbad database, is incorrect, as Ha emission was never 
reported for HR 1847 A. 



2.3. Magnitudes and photometry 

The first reli able magnitude v alues (Va - 6.07, Vb - 6.44) were 
published by CannonI (1191 2bh . She also gives the combined AB 
magnitude as 5.49. Maybe that only the combined magnitude 
value was given in her earlier spectral classification (Cannom 
Il912ah and in the HD catalogue CCannon & Pickering. .19181) 
caused its being incorrectly but often quoted as the magnitude 
of the A component until now. A similar value is also listed in 
the Bright Star Catalogue (V = 5.46, Hoffleit & Jaschek, 1991). 
The photometry of HR 1 847 was measured by several au- 
thors always during observing runs devoted to large sample of 
stars. lEggenl (1977) measured uvby/3 photometry of early type 
stars and derived values Va - 6.1 and Vb - 6.5. Lindroos 
(Il983h measured uvbyfi photometry of visual binaries includ- 
ing HR 1847. He derived corresponding V magnitudes of both 
components from the y filter values as Va - 6.091 + 0.019 
and Vb - 6.510 + 0.023. BVR observation s of visu a l bina - 
ries at Mount Maidanak were performed by IShatskiil (Il998h . 
and he obtained Va = 6.09 + 0.02 and Vg = 6.51 + 0.02. 
The EXPORT UB VRI photometry (obse rved by the Nordic 
Optical Telescope, lOudmaijer et al.L 1200 l l) gave values Va - 
6.09 + 0.10 and Vb - 6.49 ± 0.07. Recent determinations 
of magnitudes were deriv ed from HIPPARCOS photometry 
(see iPerrvman et aP. Il997l). With the HIPPARCOS Vt and Bj 
magnitudes, iKharchenko (2001) derived the magnitudes Va — 
6.056 + 0.008 and Vb = 6.451 + 0.004, which do not differ 
too much from the first values of Cannon. A more complete list 
of photometric magnitudes of both components of HR 1 847 is 
given in Table [T] 



2.6. Rotation 

Ro tational velocitie s of both components were first determined 
bv ISlettebakI (Il963h . The projected rotational velocity of com- 
ponent A was measured as (v sin /)a = 50 km s ' , while compo- 
nent B was fo und to be rotating rapidly, {u sin /)b - 300 km s ' . 
iLevatol (Il975[) determ ined (u sin /)a - 60 km s ' and (v sin Ob = 
300 km s~' . lAbt et al.1 (|2002) refined the values of the rotational 
velocities to (usini)^ = 45kms"' and (usinOB - 200km s"' 
using coude spectra obtained at the Kitt Peak 0.9-m telescope. 



2.7. Variability 

Slightly different values of magnitude obtained at different times 
using different instruments, which are listed in the Section 
12.31 cannot be considered as firm proof of variability. The 
only available homogeneous set of observations was obtained 
by the HIPPARCOS satellite. Using HIPPARCOS photometry, 
iPercv. Harlow. & Wul (12004) determined the characteristic time 
scale of variability for HR1847 as 0.9 days. They conclude 
that it is a classical Be star, but it is not clear which compo- 
nent they are refering to, since both components A and B have 
HIPPARCOS photometry available. In addition, data for both 
components also include several values that correspond to the 
combined magnitude of both components. Since it is not clear 
that they excluded the apparently wrong values (those corre- 
sponding to combined A-hB magnitude), their result is question- 
able. We performed an independent search for var iability and 
could not confirm the variability time scale found bv lPercv et al.l 
(120041) . 
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2.8. X-rays 

HR 1 847 has also been identified as an X-ray source by ROS AT 
("Bergho fer et all 1 19961) with an X-ray luminosity \ogLx - 
29.64 (in the paper version). In the online version of the same 
catalogue (their Table 2), this value is slightly different, log Lx - 
29.79. Both corresponding values of log(Lx / L\,oi) - -6.82 (or 
-6.79 in the online version) correspond to a ty pical log(Lx/Z^boi ) 
relation for B-type X-ray emitters (cf. B erghofer et al.L Il997h . 
It is, of course, a question of which star of the pair is 
the X-ray source, if not both. The position of the ROSAT 
source IRXS J053214.9-H170319 suggests it is more prob- 
ably coincident with the s tar TYC 1 301- 1942- 1=HR1847B 
dFles ch & Hardcastle'. '2004"); however, the positional error is 
18" ( Vogesetal.. 2000) , so HR 1847 A also hes within the er- 
ror circle. 



2.9. IR excess 



In their study of IR excess of 101 Be stars lCote & Waters! (Il987h 
found IRAS magnitudes [12] = 4.95, [25] = 2.22, [60] = -0.92, 
which places HR 1847 in the [12]-[25]/[12]-[60]colour-colour 
diagram far away from the region where almost all Be stars are 
located. Cote & Waters suggest ther e is a reflection nebula close 
to HR 1847. WhitelocketaLl dTgsgl) assoc iate HR 1847 with the 
IRAS SSSC source X0501-H589. Magnie retal.l(ll999h suggest 
that HR 1847 is a good candidate for a Herbig Ae/ Be star with 
a very large or a very cool circumstellar disk, but llwata et all 
(1 19991) did not find any CO emission for this star. 

Again, it is not clear which star has IR excess, but most likely 
it is component B, which was found to be a Be star There is 
also a small possibility that both stars are Be stars, with com- 
ponent A not in emis sion at that moment. F urther observations 
may shed a light on it. lCote & Watersl(ll987h attributed the spec- 
tral type B7IIIe to HR 1847, which corresponds to component B, 
and u sin / = 55 km s ', which corresponds to component A. 



2.10. Cluster membership 

Although HR 1847 is located in the same region as the open 
cluster Collinder 65 (A component is #771, B component 
#772, according to the WEBDAQ database), its membership 
is improbable mainly be cause to the different proper motion 
dKharchenko et al.Ll2004 . 



3. Observations and data reduction 

The data available for this study consist of several data sets of 
electronic spectra mostly centred on the Ha region; 



Component A 

- 68 spectra in the spectral range 6250-6770 A obtained with 
a CCD SITe ST-005 800x2000 pix camera attached to the 
coude spectrograph of the 2m-telescope in Ondfejov (Czech 
Republic). The spectra were obtained between 18 October 
2003 and 25 February 2008. Spectra were reduced using 
IRAlS 



' The Galactic and Magellanic Clouds open cluster database 
WEBDA is available at http://obswww.unige.ch/webda/ 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



- One spectrum obtained on 22 March 2003 with the red 
(5850-8450A) and blue (3800-5650A) channels of the 
fiber-fed echelle spectrograph HE ROS (re solving power ~ 
20000, for its brief description see JKaufea l 1998). The spec- 
trograph was attached to the Cassegrain focus of the 2m- 
telescope at the Ondfejov Observatory. All the basic data re- 
duction processing was done using the HEROS pipeline writ- 
ten by Otmar Stahl and Andreas Ka ufer as an extensio n of 
the basic MIDAS echelle context (see lStahl eranil995L also 
ISkoda&SlechtaLl2002h . 

- 14 spectra obtained in the spectral range (65 10-6608 A) 
at Rozhen Observatory (Bulgaria) using the coude- 
spectrograph of the 2m RCC telescope. The CCD cam- 
era Photometries AT200 with an SITe SI003AB 1024x1024 
CCD chip (24x24 jum pixel size) was used in the f/9.5 camera 
of the spectrograph to provide spectra with a spectral res- 
olution of 36000 in the Ha region. Spectra were obtained 
between 24 October 2004 and 30 August 2007 and reduced 
with the MIDAS package. 

- One spectrum in the spectral range 6290-6745 A obtained at 
O bservatorio do Pico dos Dias (OPD, Brasil) and published 
bv lVieiria et alj (|2004|) . For details and information about the 
reduction process, see the reference above. 

In addition to the elect ronic spectra, we a l so use d the radial 
velocity measurements of iPlaskett & Youna (Il919h . For com- 
parison purposes, several values of published radi al velocities 
were u sed from the works of iFrost et al. (1926) and iDuflot et al.l 
11991 . 



Component B spectra came from the same instruments as men- 
tioned in the preceding part describing spectra of component A. 

- 38 spectra obtained using the Ondfejov coude spectrograph 
between 17 October 2003 and 22 March 2008. 

- One HEROS spectrum obtained at the Ondfejov Observatory, 
one night after the HR 1847 A spectrum, on 23 March 2003. 

- Three spectra between 28 August 2007 and 30 August 2007 
obtained at the Rozhen Observatory. 

3. 1. Radial velocity measurements 

Radial velocities were measured using the code SPEFO, which 
was developed by the late Dr Jiff Horn (see also Skoda, 199q) . 
The FITS files obtained from IRAF and MIDAS were transformed 
to the SPEFO format, and then the radial velocities (RVs) were 
obtained interactively by means of the best match of the line 
profile with its mirror. The RVs were obtained for hydrogen Ha, 
He I 6678 A, Si ii 6347 A and 637 1 A fines. RV data for individ- 
ual components are described later in this paper in Sections 14.21 
andEB 

4. Component A 

4.1. The effective temperature and gravity 

The spectral region observed by HEROS covered the range 
3800-8620 A (see Figs. lA.l land lA.2b and was used to determine 
the spectroscopic effective temperature and surface gravity. The 
line spectrum of the star is quite poor. Balmer lines up to HIO 
and some infrared lines are seen. No emission is present in the 
Balmer lines and in the lines of other metal ions. 

The effective temperature and surface gravity were deter- 
mined by comparing the observed spectra to model spectra. 
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Fig. 1. The radial velocity variations for hydrogen Ha (upper 
panel) and He i 6678 A (lower panel) lines measured from the 
Ondrejov and Rozhen spectrograms as a function of HJD. The 
He I 6678 A line is only present in the Ondrejov spectra. Full 
symbols represent Ondrejov data, and open symbols stand for 
Rozhen data. All values are from Tables [ATI and lA!2l 



We used the ATLAS9 LTE line blanketed model atmospheres, 
which were calculated by Kurucz ( 1993) assuming solar metal- 
licity and microturbulent velocity 2 km s"'. All synthetic spec- 
tra were convolved (using the code R0TIN3 by I. Hubeny) with 
the Gaussian function having FWHM=0.25 to reduce the res- 
olution of the synthetic spectra to the observed resolution. We 
used a ;t'^-fitting routine to compare the whole observed spec- 
tra with the synthetic ones, and we obtained a best fit of Teff = 
(12500 ±500)Kandlog;? - 3.5 ±0.5. F ollowing the temperature 
scale of lTheodossiou & Danezisl(ll991h . these values correspond 
to spectral type B7 - B8, and from the value of the surface grav- 
ity, we estimate the luminosity class as IV - V. The quoted error 
bars correspond to the adopted steps of our model grid, which 
are 500 K in temperature and 0.5 in log g. 

Using the code R0TIN3, we calculated a grid of rotationally 
broadened spectra for 5kms ' < usini < 500 km s' with a 
step of 5 km s"'. Then we applied the x^ fitting routine again 
to determine the best rotational velocity for this star To do this, 
we u sed the spectral line Mg ii 448 1 A (recommended by iGravL 
11976 . as a line free of pressure broadening) in the fitting process. 
We determined the projected rotational velocity of HR 1847 A 
as u sin / - (40 + 3 ) km s^ ', which differs by 5 km s"' from the 
value of lAbt et alJ (l2002h . The error of our value of the rota- 
ti onal velocity was ca lculated employing a 1-cr error algorithm 
of lZhangetal](ll986h . 



Tables lAn and lA.2l displav the line radial velocity measurements 
that were derived using the methods described in Section 13.11 
Columns 1, 2, and 3 represent the file identification, HJD, and 
the heliocentric RVs, respectively. Column 4 lists the RV mea- 
surements for the Ha line. Columns 5, 6, and 7 list the RVs mea- 
surements for He I 6678 A, Si n 6347 A, and Si ii 637 1 A (which 
are weak) absorption lines, respectively. The measured RVs were 
shifted to the zero-point using a set of shar p telluric ab s orptio n 
lines by means of the technique described in lHorn et alJ (Il996l) . 

The RVs of HR 1847A vai-y between 8 km s ' to 27 km s"' 
with a mean velocity around 13km s'. Figure [T] displays the 
measured RVs for hydrogen Ha and He 1 6678 A as a function of 
time. The plots indicate that our RV measurement started with 
one value just before the i?V maximum, then after the maximum, 
a clear decrease in the RV values is recorded until the RV val- 
ues reach their minimum around HJD 2453620 with RV values 
of 9-10 km s '. Then they quickly reached the next maximum 
value of 27 km s ' around HJD 2453670. Although only two 
cycles are covered by our observations, fortunately we success- 
fully recorded two epochs of rapid rise to the maximum, which 
happened during a very short time (nearly 35 days) of the long 
cycle. Rozhen observations confirmed the present distribution of 
RVs. Such a relatively short event could easily be missed. 

The i?V measurements that were obtained with different lines 
in the Ha region show the same distribution with time, and only 
a small difference in the maximum height can be recognized, as 
illustrated for the He 1 6678 A line (lower panel of Fig.[l). 

4.3. Period determination 

One of the most fundamental parameters of binary systems is the 
orbital period. The present RVs span an interval of 1801 days. 
From the inspection of the time plot one can conclude that 

1 . the present RV measurements of Ha and of some other met- 
alic lines show clear evidence of long-term cyclic variability, 
and 

2. the relatively large scatter between individual measurements 
on successive nights indicates that shorter period variability 
could be present in the data distribution. 

No period search of the RV measurements has been carried 
out so far for HR 1847 A. Improper data distribution, gaps, and 
timing of the observations are the most common problems in the 
period search process. The period search was carried out sep- 
arately for several sets of measurements, and two independent 
numerical period searching routines were used. The first one is 
based on the phase dispersion minimization (PDM) technique 
(Stellingwerf, 1978), the other one, PERIOD dBreged [T990l) . is 
based on a Fourier analysis. 

For PERIOD, starting values for the frequencies need to be 
given, and frequency values are improved within the limits given 
by the window function. We used this program to search for 
periodicities in both Ha and Hei 6678 A RV data sets in the 



Table 2. Candidate frequencies resulting from the PDM period 
search. 



frequency (d ) period (days) 6 



0.0069 
0.0014 
0.1077 



144.92753 0.3841 

714.28571 0.4092 

9.28505 0.4668 
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Fig. 2. Upper panel: Power spectrum of measured RVs of Ho- 
and He i 6678 A in the frequency range from 0.001 d ' to Id', 
middle panel: spectral window, and lower panel: power spec- 
trum of the prewhitened residual with 0.0014d '. In all panels 
frequencies at the .x-axis are given in d '. 



range between 1 to 1000 days, with an expected frequency res- 
olution of 0.00055 d'. The upper panel of Fig. |2] displays the 
periodogram of the detected frequencies in this period range. 
It shows three candidate frequencies (indicated by arrows in 
the upper panel of Fig.|2]i, namely 0.0069d-', O.OOUd', and 
0.1077d"^ The frequency 0.0069d ' has power 61, and is 
dominant. The middle panel of Fig. |2] shows the spectral win- 
dow. In the lower panel of Fig.|2]we show the power spectrum 
of the residual periodogram prewhitened for 0.0014d"' (using 
PERIOD), which clearly shows the disappearance of all peaks re- 
lated to all three candidate frequencies. 

Using the PDM technique we searched for periodicity in in- 
tervals 1-lOdays, 1-100 days, and 1-1000 days. The results are 
summarized in Table |2] and agree with the ones obtained with 
PERIOD. Figure [3] illustrates the phase diagrams of measured 
RVs for the Ho- line folded with the 0.1077 d', 0.0069 d' and 
0.0014 d"' frequencies. 

Although the frequency 0.0069 d' has the highest power 
(cf. Fig.|2] and Table|2]l, we suspect that this frequency is not a 
real on e. First, the historical observations of iPlaskett & Young 
( Il919l) are out of phase by about half of the period (open 
squares in the middle panel of Fig. [3]). Second, the correspond- 
ing RV curve predicts a maximum of RVs near HJD 2453350, 
which was not confirmed by our observations. Values around 
~15 km s ' appear instead of maximum ones. In addition, the in- 
terval of 145 days (or approximately five months) roughly corre- 
sponds to the gap between the end of one season in March/ April 
and the next start of observations in August/September Third, 
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Fig. 3. Phase diagram of measured RVs for Ha lines folded with 
different frequencies displayed from top to bottom, 0.1077 d"' 
(9^28505), 0.00689d-' (144'?9275), and 0.0014d-' (714'?2857). 
Filled circles (•) and open circles (o) denote data obtained 
from our observations at Ondfejov and Rozhen, respectively. 
Open squares (d) and o pen triangles (7} repre sent data from 
l3<I ' 



iPlaskett & Younj (1 1 9 1 9h and lVieiria et alJ (l2004 . respectively. 



the candidate frequency 0.0069 d ' is almost exactly equal to 
5x0.0014 d"'. The relative scatter of individual i?V values folded 
with the 0.1077d"' frequency is quite large. Following these ar- 
guments, we propose the frequency at 0.0014d"' as the domi- 
nant frequency. 

4.4. Orbital soMion 

The variability of the radial velocities with frequency 0.0014d ' 
can most li kely be explained by binar ity, and it supports the ear- 
lier idea of IPlaskett & Yound (Il919l) that the star is a spectro- 
scopic binary. 

Starting from the frequency 0.0014 d"', we performed 4 dif- 
ferent or bital solutions using the program FOTEL developed by 
iHadraval (1990, 2004). For the solutions denoted as I, II, and 
III (see Table [3), the orbital elements were derived for Wa ra- 



J. Kubat et al.: Spectroscopic analysis of the B/Be visual binary HR 1847 



05 

E 

> 

cc 



25 












^ 


15 
5 


- 1 


■| s 




i 

- , 


•is . ' - ' 


■| S . ' 


-5 




: 


h' 


= ■1 


:'h''^ 


:.\v=i 


15 












- 



a -5 



-0.5 



0.5 



1 "I * • . ! oS 



■J.' 



-0.5 



0.5 
Phase 



1.5 



. J." 



1.5 



Fig. 4. Upper panel: Radial velocity curve of the primary 
component corresponding to the FOTEL orbital solution III 
{P = 719^79), which was based on Ho- RV measurements of 
Ondfejov, Rozhen, and OPD spectra and on RV values from 
[Plaskett & Young ( 1919). i?V values are displayed in the centre- 
of-mass reference system. Lower panel: Phase plot of O - C 
deviation from the RV curve for Solution III. 
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Fig. 5. The evolution of the Ra profile of HR 1847B. 



dial velocity variations with data observed from different sites. 
Solution I is only performed with Ondfejov data. Solution II 
is obtained using Ha measurements from both Ondfejov and 
Rozhen observations. Solution III is obtained using all of the 
available data, namely Ondfejov and Rozhen observations, one 
spectrum from OPD, and including the historical measurements. 
Solution IV was obtained using radial velocities of all available 
spectral lines in the Ha region (Ho-, He i 6678 A, Si ii 6347 A, 
and Si II 6371 A). 

In our calculations we allowed the period, eccentricity, peri- 
astron longitude, and semiamplitude to converge. For each data 
set FOTEL also allows individual y-velocities to be determined. 
Figure |4]illustrates the FOTEL Solution III with its (9 - C residu- 
als. We searched for additional periodicities in the residuals but 
did not find any significant periods. 

We adopt the Solution III (P = 719^79+(>?17,e = 0.70+0.02) 
as an orbital solution of HR 1847 A. 



5. Component B 

The B component is a Be star with a relatively weak emission 
in Ho- and almost negligible one in H/3. Since 2003, the weak 
emission in Ho- further weakened, indicating an approaching end 
of the Be phase. No time scale may be given for this appar- 
ently long-term variability, since only 4 years of observations 
are available. Long term variability (strengthening and weaken- 
ing of emission) is typical of Be s tars on t i me sc ales from years 
to tens of years (for a review see iHubertl l2007h . For example, 
X Dra has quite a well-established long term ti me scale, and 
its m ost recent determination is 22.11 years (see iKubat et al.L 
12010 . and references therein). The Ha profiles obtained by our 
observations are plotted in Figure |5] We fitted the spectrum 
of component B using the same methods as for component A. 



500) K, logg = (3.5 +0.5), and 



We found T^g = (12500 
ysin/ = (230 + 5) km s"'. 

5. 1 . Radial velocities 



Owing to the high rotational broadening of lines in the com- 
ponent B spectrum, the helium and silicon lines, which were 
used for RV measurement of component A, are too shallow 
in component B, and therefore unusable for RV measurement. 
Consequently, we measured the RVs of component B only from 
the Ha line. RVs were measured with the profile inversion tech- 
nique. However, the rapid rotation of the component B and rela- 
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Fig. 6. Time evolution of the Ha radial velocities of HR 1847B. 
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lively low S/N made it impossibile to arrive at precise RV values. 
Cross-correlation techniques, such as least-square deconvolution 
(LSD, Iponati et al., 1997) would not improve the RV accuracy 
as too few line profiles are available in our spectra. Also, the bi- 
sector method does not give more accurate results because of the 
limited quality and not high enough S/N of our spectra. 

The Ha line in the B component is characterized by a double- 
peak emission with peaks of equal strength. We measured the 
RVs of three different features at the Ha line, namely of the line 
wings, of the emission part, and of the central absorption. The 
results are listed in Tables lA3] and lA!4l and plotted in Fig.|6] 

Although the RVs are evidently variable, a period search did 
not give any reasonable result. That is why we cannot make 
any conclusions about the possible binarity of HR 1847B until 
a more complete set of more accurate data has been obtained. 



6. The possible A+B system 

There is a question as to whether the components of this visual 
pair are physically bound. However, there is not enough data to 
help to resolve this question. The history of all relative position 
measurements of HR 1847 in the sky (see Section l2.1.1b do not 
allow any conclusion about their relative motion. Slightly differ- 
ent values can be attributed to errors in individual measurements. 

Unfortunately, only a common parallax measurement from 
HIPPARCOS is available. Since individual measurements are 
missing, we cannot easily say whether the A and B components 
are at the same distance. However, study of visual binaries by 
IShatskiil (11998 !) suggests that HR 1847 could be a physical sys- 
tem (with a period shorter than IMyr), since both components 
have similar proper motion and corresponding photometric par- 
allaxes, which agree with the hypothetical parallax. 

In case the visual components of HR1847 are physically 
bound, their similar spectral type suggests that they are of sim- 
ilar age. Then there is the question of why one of these stars 
developed an envelope that causes Ha emission and the other 
does not. The emission could not be caused by mutual interac- 
tion between components A and B, simply because the distance 
between the components is too large. A different mechanism for 
the origin of the Be ph enomenon in HR 1847B h as to be found 
(see, e.g., the review bv lPorter & Riviniusll2003l) . The emission 
in the component B is most probably connected with its rapid 
rotation. Projected rotational velocity of component A is lower, 
and this star also has no emission. On the other hand, since the 
inclination angles are unknown, we can not exclude the possibil- 
ity that the component A is also rapidly rotating, as there is ap- 
parently no correlation between directions of spin axes for wide 
binaries (.Howe & Clarke. ,2009). 



6.1. Comparison with other visual binaries 

HR 1 847 is not the only example of a visual multiple system with 
similar components, with some of them also showing the Be phe- 
nomenon. Thejdsu2L£iEl5_5X5l^™_^^on consists of three Be 
stars (see iMaranon di Leo et al.L 11994 and references therein), 
and all these three stars are rapidly rotating (jAbt et al., 2002). 
The projected rotational velocities of /3 Mon A and /3 Mon C 
was reported as 260 km s"' and 250 km s"', respectively, while 
y6 Mon B is rotating more sl owly (140km s""') and, in terestingly, 
also has weakest emission ( Cowlev & GugulaLll973 ). This sup- 
ports the connection between rapid rotation and the Be phe- 
nomeno n. More visual binarie s with Be-components may be 
found in lAbt & Cardonal (11984 . 



Table 4. Summary of parameters of both components of 

HR1847. 



Q'(J2000) 
(5(J2000) 

V 



logg 

vsin; (km s"') 

Be star 

binary 



HR1847A 



HR1847B 



05 3214.14 

-H17 03 29.3 

6.09 



05 32 14.56 

+ 170321.8 

6.51 



12500 ± 500 12500 + 500 

3.5 ±0.5 3.5 ±0.5 

40 ± 3 230 ± 5 



no 
SBl* 



yes 
SBl? 



Binary parameters are listed in Table [3J 



Investigation of such multiple systems with resolved visual 
components (like HR 1847 or /3 Mon) may help when study- 
ing multiple partially resolved or unresolved systems with Be 
stars like o And (see Olevic & Cvetkovic, 2006, and references 
therein) or /3 Cep (see Whe elwright et al.L 1.2009, and references 
therein), where light from all stars is mixed, making the analysis 
extremely difficult. The physical parameters and system charac- 
teristics derived from resolved systems will put constraints on 
stellar structure and binary evolution models, such that the un- 
resolved systems can be modelled in a better way with the im- 
proved models. Careful study of resolved systems and careful 
determination of their parameters are very important feedback 
on models. 



7. Conclusions 

This paper is the first attempt to collect all available information 
about the stars in the visual binary system HR 1847 and to cor- 
rect errors that appear in the SIMBAD database and Bright Star 
Catalogue. We have present ed results of ou r spectroscopic anal- 
ysis. HR 1847A (V = 6.09. IShatskiil 11998 ) is a B type single- 
lined, eccentric, (e = 0.70 + 0.02), spectroscopic binary with 
a period of 719.79 + 0.17 days. Model atmosphere analysis of 
this star yielded T^ff = 12500 K, log.? - 3.5, an d (ysinOA = 
40 km s'. HR1847B (V = 6.51. IShatskiil [19981) is a Be star 
with variable radial velocities, however, no reasonable orbital 
solution could be found based on the available data. Model at- 
mosphere analysis of this star yielded the same effective temper- 
ature and surface gravity as for the A component, but with a ro- 
tational velocity of (v sin Ob - 230 km s"' . A decrease in the Ha 
emission was recorded during the last 4.5 years of observations. 
In addition, positions of X-ray (IRXS J053214.9+170319, see 
Section l278b and infrared (X0501+589, see Section l2!9] l sources 
indicate that HR 1847B is more likely an X-ray and IR source 
than HR 1847 A. 

Several other issues, such as the short term variability of 
HR1847A or the reliable period determination of HR1847B, 
could not be addressed by this paper owing to the lack of neces- 
sary data. Consequently, further long-term observations of these 
stars are desirable. Future work should also include a detailed 
abundance analysis using high S/N spectra and NLTE model at- 
mospheres. 
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Table 3. Orbital elements of HR 1847A with FOTEL. 



Element 



Solution I 



Solution II 



Solution III 



Solution IV 



P[d] 


715.03 ±1.06 


713.83 ±1.49 


719.79 ±0.17 


710.49 ±0.98 


^ periast. 


51501.61 ±4.26 


51504.34 ±5.93 


51480.31 ±3.88 


51503.67 ±3.09 


A:(kms-') 


7.89 ±0.53 


7.58 ±0.55 


7.73 ± 0.44 


8.23 ± 0.23 


e 


0.79 ± 0.02 


0.72 ± 0.03 


0.70 ± 0.02 


0.73 ± 0.01 



w[deg] 



306 ± 5.87 



310 ±6.99 



304 ± 7.03 



286 ± 3.70 



yi(km s'^ 
-y2(krn s" 
73(1011 s" 
y4(km s" 
ysCkm s'^ 
yeCkms" 
77(km s" 
ysCkms" 

f{m) 



13.47 



0.0095 



3.39 


13.44 


- 


9.54 


9.88 


- 


- 


10.63 


— 


- 


10.52 


- 


- 


- 


13.23 


- 


— 


13.36 


- 


— 


10.97 


_ 


_ 


13.10 



0.0104 



0.0126 



0.0132 



No. ofRVs 


64 


76 


83 


246 


rms (km s"') 


1.76 


2.08 


2.22 


1.94 



Note: The velocity 71 is based on Ondfejov data, 72 on Rozhen data, 73 on OPD data jVieiria et all 120041) . and 74 on historical data of 
iPlaskett & Yound dl919l) . In solution IV 75, 75, 77, and yg stand for Ha, He i 6678 A, Si 11 6347 A, and Si 11 6371 A, respectively. 
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Appendix A: Visual spectra of IHR 1847A and 
HR1847B. 



Table A.2. Heliocentric radial velocities (RV) of HR 1847 A 
spectra from Rozhen Observatory. 



file 


HJD 


Vhel 


RV(Ha) 




(JD-2400000) 


(kms"') 


(kms-') 


04 j 119 


53303.4455 


23.17 


10.52 


04J12O 


53303.4563 


23.16 


9.14 


04 j 303 


53306.5530 


21.99 


11.16 


04 j 304 


53306.5638 


21.97 


12.06 


04kl09 


53333.3218 


10.99 


8.39 


05C062 


53452.6889 


-29.64 


4.76 


05cl63 


53453.3861 


-29.56 


7.58 


05c235 


53454.4102 


-29.47 


4.93 


07h858 


54341.5910 


28.25 


11.92 


07hl29 


54342.5379 


28.43 


12.12 


07hl35 


54342.6011 


28.36 


12.05 


07h215 


54343.5844 


28.52 


12.19 



relative intensity 



relative intensity 



relative intensity 



relative intensity 
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Table A.l. Heliocentric radial velocities (RV) of HR 1847A - spectra from Ondfejov Observatory, where Vhei is the heliocentric 
velocity. 



file 


HJD 


Vhd 


RV(Ha) 


RV (He I 6678 A) 


RV (Si II 6347 A) 


RV(Siii6371A) 




(JD-2400000) 


(kms-') 


(Ions-') 


(kms-') 


(kms-') 


(kms-') 


ra5788 


52721.3728 


-29.71 










mj 170023 


52930.5671 


25.36 


21.45 


21.11 


18.65 


20.07 


mj 170025 


52930.6315 


25.23 


21.03 


20.75 


17.80 


18.98 


mkl20O38 


52956.5608 


16.02 


21.45 


21.67 


17.67 


20.76 


ml090O33 


52983.4198 


3.03 


20.12 


18.98 


16.15 


18.86 


na050O19 


53010.4552 


-11.13 


13.55 


15.33 


14.84 


16.80 


na050O21 


53010.5053 


-11.23 


16.48 


15.82 


13.41 


16.43 


nd210O18 


53117.2973 


-23.56 


13.44 


16.89 


15.73 


16.70 


nd210O19 


53117.3210 


-23.57 


15.78 


17.09 


13.39 


12.61 


nd220O14 


53118.3330 


-23.35 


15.06 


18.36 


13.79 


16.32 


nd3O0O02 


53126.3175 


-20.47 


12.67 


13.76 


14.90 


18.87 


ng3O0O48 


53217.5797 


21.00 


18.70 


10.85 


15.09 


18.95 


ng3O0O49 


53217.5864 


21.01 


16.35 


12.35 


12.79 


17.14 


nh220O44 


53240.5571 


27.37 


15.60 


13.35 


14.20 


13.74 


nh220O45 


53240.5649 


27.36 


14.43 


13.58 


13.35 


15.67 


nh220O49 


53240.5862 


27.35 


15.00 


14.95 


12.49 


14.69 


nh220O53 


53240.6052 


26.32 


14.99 


10.90 


13.32 


15.52 


nh230O63 


53241.5625 


27.55 


12.98 


14.95 


13.63 


17.10 


ni020O57 


53251.6021 


28.97 


15.21 


15.56 


11.75 


15.84 


ni060O34 


53255.6179 


28.14 


11.90 


11.14 


9.90 


14.66 


ni080O43 


53257.5953 


29.47 


14.39 


— 


14.39 


— 


nj 240062 


53303.5864 


23.00 


11.76 


10.54 


9.57 


12.05 


nk240O49 


53334.5637 


10.01 


8.79 


— 


— 


— 


nllO0O2O 


53350.2773 


2.36 


13.55 


11.30 


9.52 


11.50 


oal60O18 


53387.3315 


-16.40 


12.44 


12.31 


10.18 


11.89 


OC220O1O 


53452.3157 


-29.61 


10.35 


10.01 


6.21 


5.04 


OC310O22 


53461.3659 


-28.62 


10.47 


10.03 


6.91 


10.42 


od010O15 


53462.2998 


-28.42 


11.33 


9.78 


6.25 


11.64 


od010O16 


53462.3100 


-28.41 


10.15 


10.12 


5.51 


10.18 


01060036 


53620.5833 


29.34 


9.90 


9.15 


7.15 


4.46 


01060037 


53620.5954 


29.33 


9.30 


11.79 


8.23 


10.97 


oj 0700 18 


53651.4465 


27.84 


23.89 


21.65 


20.52 


20.50 


OJ070O28 


53651.4953 


27.79 


23.28 


21.63 


21.83 


23.74 


OJ070O34 


53651.5630 


27.69 


24.35 


21.76 


20.64 


22.79 


oj 070040 


53651.6182 


27.58 


23.69 


21.45 


19.95 


21.99 


oj 070042 


53651.6418 


27.54 


23.05 


23.36 


20.76 


21.94 


oj 070044 


53651.6703 


27.48 


25.93 


24.46 


21.79 


23.70 


oj 0800 5 8 


53652.4732 


27.62 


25.34 


26.92 


25.08 


25.40 


oj 080062 


53652.5220 


27.56 


23.91 


23.42 


21.44 


23.02 


oj 080069 


53652.5815 


27.46 


21.98 


23.37 


20.87 


22.44 


OJ080O73 


53652.6229 


27.38 


24.64 


28.68 


25.40 


23.28 


OJ080O79 


53652.6823 


27.26 


27.72 


25.32 


23.06 


27.61 


oj 260033 


53670.4313 


21.60 


26.45 


24.66 


19.67 


23.76 


oj 26003 5 


53670.4447 


21.58 


24.58 


22.92 


21.83 


22.17 


oj 290036 


53673.4859 


20.51 


25.52 


22.75 


20.12 


21.90 


oj 290037 


53673.5018 


20.48 


21.98 


22.49 


20.33 


22.87 


pa080O32 


53744.3096 


-12.56 


16.11 


16.95 


14.12 


14.79 


pj 110024 


54020.6019 


26.67 


11.60 


11.71 


9.29 


12.49 


qcl60O19 


54176.2664 


-29.87 


9.03 


10.00 


5.16 


8.53 


qcl60O37 


54176.3253 


-29.98 


8.83 


10.50 


6.73 


7.65 


qcl60O46 


54176.4145 


-30.04 


9.94 


9.52 


6.06 


8.80 


qc270Oll 


54187.3995 


-29.25 


8.95 


10.44 


5.00 


7.82 


qd020OlO 


54193.3620 


-28.37 


7.03 


10.44 


4.54 


7.95 


qd070O02 


54198.3054 


-27.40 


9.50 


10.71 


6.31 


8.86 


qdl30O19 


54204.2939 


-26.03 


8.63 


9.98 


3.03 


7.47 


qdl30O23 


54204.3337 


-26.05 


9.20 


9.15 


5.43 


4.80 


qdl40O22 


54205.3262 


-25.79 


9.31 


9.81 


5.81 


8.30 


qe020Oll 


54223.3037 


-20.03 


13.92 


— 


— 


— 


qe020O12 


54223.3097 


-20.03 


11.02 


12.46 


2.89 


6.97 


qh240O29 


54337.5997 


27.59 


11.53 


13.34 


12.39 


15.14 


qh2 50033 


54338.6177 


27.71 


11.84 


11.91 


9.87 


13.95 


qll50O36 


54359.5365 


29.70 


20.28 


22.21 


18.66 


20.35 


qll60O31 


54360.5360 


29.70 


20.50 


22.06 


19.61 


23.03 


qll60O39 


54360.6371 


29.57 


21.48 


20.63 


17.17 


21.79 


qj 130036 


54387.4787 


26.51 


21.26 


18.68 


21.17 


20.95 


qj 130040 


54387.5534 


26.39 


22.08 


20.12 


20.21 


21.70 


qj 130046 


54387.6826 


26.13 


23.03 


21.13 


19.13 


21.96 


rb250O29 


54522.3202 


-28.60 


17.53 


16.29 


15.33 


17.04 
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Table A.3. Heliocentric radial velocities (RV) of HR 1847B - spectra from Ondfejov Observatory. 



file 


HJD 


Vhel 




RV(lia) 






(JD-2400000) 


(kms-') 


(kms-') 


(kms-') 


(kms-') 








wings 


emission 


absorption 


ra5806 


52721.8197 


-29.59 


11.10 


9.73 


14.30 


mj 17802 5 


52930.6315 


25.23 


15.40 


12.47 


18.92 


mkl2Q038 


52956.5608 


16.02 


15.05 


12.71 


19.70 


nc290015 


53094.3765 


-28.89 


14.39 


13.21 


17.90 


ndSOQOOl 


53126.3043 


-20.50 


18.17 


16.41 


18.17 


nh220O47 


53240.5755 


27.33 


27.44 


16.00 


12.75 


ni020O58 


53251.6217 


27.36 


7.47 


6.88 


9.22 


nil30836 


53262.5828 


29.65 


10.24 


9.06 


10.82 


nj 240860 


53303.5565 


23.02 


3.83 


-0.19 


2.16 


nj 240861 


53303.5692 


22.99 


8.58 


7.40 


7.99 


oal60816 


53387.3113 


-16.39 


19.80 


8.08 


18.92 


od010817 


53462.3213 


-28.44 


10.92 


9.16 


10.92 


01060838 


53620.6195 


29.30 


21.00 


18.06 


24.51 


01220834 


53636.6211 


29.34 


23.59 


24.17 


24.76 


01220835 


53636.6365 


29.31 


17.11 


20.04 


23.52 


OJ070834 


53651.5630 


27.69 


19.70 


17.94 


29.08 


oj 070848 


53651.6182 


27.58 


12.55 


-2.11 


13.14 


oj 070842 


53651.6418 


27.54 


21.86 


30.69 


21.89 


oj 070844 


53651.6703 


27.48 


21.25 


14.80 


27.70 


oj 080868 


53652.4977 


27.60 


— 


37.52 


24.62 


oj 280842 


53672.5077 


21.87 


12.90 


16.42 


19.94 


oj 29083 5 


53673.4623 


21.60 


11.50 


7.39 


— 


pj 110826 


54020.6245 


26.62 


7.44 


19.17 


13.31 


qcl60819 


54176.2664 


-29.87 


8.43 


7.26 


10.78 


qcl60839 


54176.3427 


-29.98 


7.94 


5.60 


11.46 


qcl60848 


54176.4353 


-30.05 


8.46 


7.87 


11.39 


qc270818 


54187.3763 


-29.24 


2.51 


7.78 


11.30 


qd020811 


54193.3770 


-28.37 


5.14 


22.14 


3.97 


qdl30821 


54204.3139 


-26.05 


8.04 


3.93 


11.56 


qdl40824 


54205.3461 


-25.80 


4.30 


8.99 


7.23 


qh240838 


54337.6245 


27.57 


20.89 


32.61 


23.23 


qil50838 


54359.5627 


29.68 


22.60 


22.02 


21.43 


qil60833 


54360.5562 


29.68 


23.68 


37.17 


24.27 


qil60839 


54360.6371 


29.57 


20.12 


37.13 


20.12 


qj 130846 


54387.6826 


26.13 


22.66 


32.04 


23.25 


rb250827 


54522.2805 


-28.54 


8.89 






rb250833 


54522.3822 


-28.73 


8.21 






rb250835 


54522.4158 


-28.77 


7.35 






rc220812 


54548.3622 


-29.65 


21.58 







Table A.4. Heliocentric radial velocities (RV) of HR 1847B - spectra from Rozhen Observatory. 



file 


HJD 


Vhel 




J?V(Hq') 






(JD-2400000) 


(kms-') 


(kms-') 


(kms-') 


(kms-') 








wings 


emission 


absorption 


07h849 


54341.5732 


28.25 


21.98 


35.69 


23.52 


07hl36 


54342.6176 


28.33 


21.19 


28.47 


26.19 


07h216 


54343.5993 


28.50 


21.33 


25.90 


21.33 



